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Abstract 

We discuss a formulation of exactly Poincare invariant quantum me- 
chanics where the input is model Euclidean Green functions or their gener- 
ating functional. We discuss the structure of the models, the construction 
of the Hilbert space, the construction and transformation properties of 
single-particle states, and the construction of GeV scale transition matrix 
elements. A simple model is utilized to demonstrate the feasibility of this 
approach. 



1 Introduction 

The motivation for this work is to construct mathematically well-defined quan- 
tum mechanical models of few-body systems at the GeV energy scale that have 
a direct relation to an underlying quantum field theory. The goal is to use ex- 
perience gained from the field theory to constrain the structure of the models. 

We do this by starting with the quantum mechanical interpretation of the 
field theory. Normally this is given in terms of vacuum expectation values 
of products of fields (Wightman functions), which represent the kernel of the 
Hilbert space inner product of the field theory [5] ; however the Wightman func- 
tions do not have a simple connection with the Lagrangian formulation of the 
field theory. The Euclidean Green functions have the advantage that they can 
be directly related to Lagrangian field theory and at the same time can be used 
to reconstruct the underlying quantum theory [3] @] [5] . 

With this connection in mind we consider a class of models that are ideally 
expressed in terms of Euclidean-invariant reflection-positive Green functions or 
their generating functional. The generating functional is formally the functional 
Fourier transform of the path measure: 

Z[f] := J D e [ct>]e- m+iW) = £ ^G n (fr^J) = ]J exp gc„(/, •••,/) 

n times 

(1) 
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This provides the formal relation to the field-theoretic Lagrangian. For the pur- 
pose of illustration we restrict our considerations to generating functionals for 
scalar fields. The G n are the n-point Euclidean Green functions smeared over 
test functions in Euclidean space-time variables and the C n are the correspond- 
ing connected Green functions. 

The generating functionals are assumed to be Euclidean invariant, reflec- 
tion positive, and satisfy space-like cluster properties. The test functions are 
assumed to be Schwartz functions with support for positive Euclidean time. We 
denote this space by S + . Euclidean transforms on the test functions are denoted 
by /(x) — > /o,a(x) := /(Ox + a)) and Euclidean time reflection is denoted by 
9/(t, x) := /(— r, x). The requirements on the generating functional are JSJ[5]: 

Z\f\ = Z[fo. a ] Euclidean invariance (2) 

6 S + Mij = Z[fi - Qfj) > reflection positivity (3) 
lim (Z[f + gi iBL ] - Z[f]Z[g}) ->• cluster properties. (4) 

a|— >oo 

Models can be constructed by specifying a set of model connected n-point func- 
tions, C n in (JTJ). A sufficient condition for reflection positivity is that each term 
in the product ([1]) is reflection positive. 

A dense set of vectors in the model Hilbert space is given by functionals of 
the form: 

N b N c 

B[4>] = 2 ftj-e**^ C[4>] = J2 Cfce 80(9fc) (5) 
j=i k=i 

where bj, Ck £ C, fj, gu £ S+ and A^, N c < oo. The model Hilbert inner product 
of two such vectors is 

(B\C) : Y^KnZg, H./V = / D e [<l>]e- A MB*[4>oe]C[<j>]. (6) 

3=1 fe=l J 

The representation at the end of eq. © suggests that we can think of the 
vectors as wave functionals, however the computation of the inner product only 
requires the generating functional. The reflection positivity condition ensures 
that vectors have positive length: 

(B\B) > 0. (7) 

To understand how Poincare invariance is realized observe that the determi- 
nants of the following matrices are (— ) the squares of the Lorentz and Euclidean 
lengths respectively: 

x _ ( t + z x-iy \ x = f ir + z x-iy \ , g . 
\ x + iy t — z J \ x + iy it — z J 

The group of linear transformation that preserves both of these determinants is 
ST(2,C) x ST(2,C): 

X' = AXB t X = AXB l dct(A) = dct(S) = 1. (9) 
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These represent complex Lorentz or complex orthogonal transformations. Real 
Lorentz transformations have B = A* while real orthogonal transformations 
have A,BG SU{2). The group of real orthogonal transformations form a sub- 
group of the complex Lorentz group in the inner product ([5]) . When one accounts 
for the support condition on the test functions, Euclidean time evolution be- 
comes a contractive semigroup, rotations in Euclidean space-time planes become 
local symmetric semigroups [7] [8] [9], and translations in a fixed direction and 
rotations about a fixed axis become unitary one-parameter groups. The gener- 
ators of these transformations are self-adjoint operators on the physical Hilbert 
space. The one-parameter groups (semigroups) can be expressed in terms of 
their infinitesimal generators as 

e- f3H f3>0 e iaP e iJ '^ e*^' . (10) 

It is straightforward to show that the generators {H, P, J,K} satisfy the com- 
mutation relations of the Poincarc Lie algebra. 

In this framework particles are point spectrum eigenstates of the square of 
the mass operator: M 2 := H 2 — P 2 . Normalizable mass eigenstates can be 
represented as wave functionals 

B,M = ^6 n eW"l (11) 

n 

Simultaneous eigenfunctionals of mass, linear momentum and angular momen- 
tum can be extracted from these mass eigenstates using rotations and transla- 
tions: 

BxfpM = J j^e-^B XJ ,M (12) 
Ba, 3 (p,m)M := / dR V B^oiR-'pMDfriR). (13) 

JSU(2) „t^j 

The single-particle wave functionals can be interpreted as multiplication op- 
erators. These single-particle wave functionals can be used to construct the 
two Hilbert space injection operators that define the asymptotic conditions in 
the two Hilbert space [TU] formulation of Haag-Ruelle Scattering theory [TT] [12"] 
|13] | 14 j . The wave operators and injection operator have the form 

|*±(/i, ' ' ' /»)) := I™ e^e^'lf) = fi±|f) (14) 

( — > DO 

$|f>M = 

/ ^n^kiPkjBxkjkiPk, fJ>k)[4>] - [H, B XkJk (p k , VkM}) fk(Pk, Vk)dp k - 
J k 

(15) 

The asymptotic Hilbert space is the tensor product of one-particle irreducible 
representation spaces of the Poincare group. Existence of the wave operators 
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can be checked in a given model by verifying the finiteness of the integral [15] : 



/>±oo 

/ \\(H<I> - <S>H )e- mot \f)\\dt < oo, (16) 

while Poincarc covariancc of the wave operators, 

[/(A, o)Ot = O±J7 (A, a) (17) 

can be checked by verifying the asymptotic condition in this representation of 
the Hilbcrt space 

lim ||(K$-$K )e- lHot |f)|| = 0. (18) 

Approximate sharp-momentum transition matrix elements can be computed 
from S matrix elements in narrow wave packets using |16) 

/ / / i i \rn\ \ (^f\S\fi) — 5 a b{ff\fi) 



Using the Kato-Birman invariance principle [17] [18] [13] [14] the expression for 
the wave operators can be replaced by the limits 

0±:= lim e - lHt <f>e lHot = lim ^ e ~^ m \ (20) 

t— >±oo n— >±oo 

which for large enough n gives the approximate expression for the S'-matrix 
elements in normalizable states: 

(f f \S\fi) = (f/|Q^-|fi) « (f/le-^^^^te^^^^e-™ 6 ""^!^). (21) 

The compactness of the spectrum of cxp(— j3H) means that for large but fixed 
n that e 2me can be uniformly approximated by polynomial in exp(— f3H): 



g 2me~ /3 - ff 



J2cm(n)(e-P mH ). (22) 



Chebyshcv expansions provide an accurate approximation |19| for large values 
of n: 



f(e-n » \c,T Q {e-n + J2 c k T k (e~n (23) 



k=l 

2 A „, ,2fc-l7i\ ,.2fc-l7i\ 



/cos )cos(7 . (24 

J \ y N+ i2> N + 1 2 V y 

We demonstrate the feasibiliy of this computational method using an exactly 
solvable relativistic model with a mass square operator given by 

M 2 = 4(k 2 +m 2 )-|. 9 )A( 5 | (25) 
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^ = <Tk* (26) 

where m is mass of a nucleon and A is chosen to give a bound state with the 
mass of a dcuteron. First we test the approximation in equation (|19p . We 
use Gaussian wave packets of the form e - a ( k -ko) and find that to get sharp- 
momentum transition matrix elements to a 0.1% accuracy the width of the wave 
packet, k w = needs to be about 3% of the initial momentum, fco- This 

works at least up to 2 GeV. The results are illustrated in table 1: 

Table 1 



fco 


a 


h 


% error 




[GeV] 


[GeV- 2 ] 


[GeV] 






0.1 


105000 


0.00308607 


0.1 


0.030 


0.3 


10500 


0.009759 


0.1 


0.032 


0.5 


3000 


0.0182574 


0.1 


0.036 


0.7 


1350 


0.0272166 


0.1 


0.038 


0.9 


750 


0.0365148 


0.1 


0.040 


1.1 


475 


0.0458831 


0.1 


0.041 


1.3 


330 


0.0550482 


0.1 


0.042 


1.5 


250 


0.0632456 


0.1 


0.042 


1.7 


190 


0.0725476 


0.1 


0.042 


1.9 


150 


0.0816497 


0.1 


0.042 



Next we test the approximation in (|21[) for the wave packet widths in table 1. We 
choose j3 so j3 times the center of momentum (CM) energy is a number of order 
unity. Table 2 shows that for n = 300 we get ten figure accuracy in the real and 
imaginary parts of the S'-matrix elements for a 2GeV incident CM momentum. 
Similar results are obtained for all momentum scales between 100 MeV and 1.9 

Table 2: fc = 2.0[GeV], a = 135[GeV~ 2 ] 



n Re (0| (S n -I) |0) 



ox 



Im (0|(5 n -/)|0) 



50 -2.60094316473225e-6 

100 -2.82916859895010c-6 

150 -2.83171624670953e-6 

200 -2.83165946257657c-6 

250 -2.83165905312632e-6 

300 -2.83165905257121O-6 

350 -2.83165905190508e-6 

400 -2.83165905234917e-6 



1.94120750171791e-3 
2.35553585404449e-3 
2.37471383801820e-3 
2.37492460997990e-3 
2.37492527186858e-3 
2.37492527262432e-3 
2.37492527262493e-3 
2.37492527262540e-3 



-2.83165905227843c-6 2.37492527259701e-3 



GeV. 

Finally we test the Cheybshev approximation for the wave packet widths in ta- 
ble 1 and the n- values in table 2. Table 3 shows that for polynomials of degree 
slightly larger than n one obtains a 10-13 figure accuracy uniformly for spectrum 
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of exp(-pH). 



Table 3: Convergence with respect to Polynomial degree e 



X 


n 


dcg 


poly error % 


0.1 


200 


200 


3.276e+00 


0.1 


200 


250 


1.925c-ll 


0.1 


200 


300 


4.903e-13 


0.1 


630 


630 


2.069e+00 


0.1 


630 


680 


5.015e-08 


0.1 


630 


700 


7.456e-ll 


0.5 


200 


200 


1.627e-13 


0.5 


200 


250 


3.266e-13 


0.5 


630 


580 


1.430e-14 


0.5 


630 


680 


9.330e-13 


0.9 


200 


200 


3.276e+00 


0.9 


200 


250 


1.950e-ll 


0.9 


200 


300 


9.828e-13 


0.9 


630 


630 


2.069c+00 


0.9 


630 


680 


5.015c-08 


0.9 


630 


700 


7.230c-ll 



Table 4 shows the final approximation for the real and imaginary parts of the 
sharp- momentum transition matrix elements for CM momenta up to 1.9 GeV. 
The results are all within less than 0.1% of the exact results in this model. 

Table 4: Approximate transition matrix elements 



k 


Real T 


Im T 


% error 


0.1 


-2.30337e-l 


-4.09325e-l 


0.0956 


0.3 


-3.46973e-2 


-6.97209e-3 


0.0966 


0.5 


-6.44255e-3 


-3.86459e-4 


0.0986 


0.7 


-1.88847c-3 


-4.63489e-5 


0.0977 


0.9 


-7.28609e-4 


-8.86653e-6 


0.0982 


1.1 


-3.35731c-4 


-2.30067e-6 


0.0987 


1.3 


-1.74947e-4 


-7.38285e-7 


0.0985 


1.5 


-9.97346e-5 


-2.76849e-7 


0.0956 


1.7 


-6.08794e-5 


-1.16909e-7 


0.0964 


1.9 


-3.92110e-5 


-5.42037e-8 


0.0967 



2 Conclusion 

We presented a formulation of relativistic quantum mechanics [1] that uses Eu- 
clidean generating functionals or Green functions as input. In applications these 
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have to be modeled. One virtue of this representation is that all calculations 
can be performed without analytic continuation. 

The expression in equation (TTJ) suggests that the generating functionals can 
be modeled using a finite collection of model connected Green functions. While 
it is easy to maintain Euclidean covariance and cluster properties of the models 
in this representation, reflection positivity is a non-trivial constraint that will 
be the subject of future investigations. While it holds for free field generating 
functionals, it is not stable with respect to small perturbations [20] . Failure of 
reflection positivity points to violations of the spectral condition or the positivity 
of the Hilbert space norms. 

The model calculations presented suggest that for models based on reflec- 
tion positive generating functionals this framework can be used to accurately 
compute both bound state and scattering observables. 

This work supported in part by the U.S. Department of Energy, under con- 
tract DE-FG02-86ER40286. 
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